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Abstract 

Motivated by the solar composition problem and by using the recently developed Linear Solar Model 
approach [I] , we analyze the role of opacity and metals in the sun. After a brief discussion of the relation 
between the effects produced by a variation of composition and those produced by a modification of the 
radiative opacity, we calculate numerically the opacity kernels that, in a linear approximation, relate an 
arbitrary opacity variation to the corresponding modification of the solar observable properties. We use these 
opacity kernels to discuss the present constraints on opacity (and composition) provided by helioseismic and 
solar neutrino data. 



1 Introduction 



In the last few years a new solar problem has emerged. Recent determinations of the photospheric heavy 
element abundances [21 [3l |3| indicate that the sun metallicity is lower than previously assumed [5j . Solar 
models that incorporate these lower abundances are no more able to reproduce the helioseismic results. 
As an example, the sound speed predicted by standard solar models (SSMs) implementing the heavy 
element admixture of [2] disagrees at the bottom of the convective envelope by ~ 10 cr with the value 
inferred by helioseismic data (see e.g. |6] and black dashed line in fig{T]). In addition, the predicted 
surface helium abundance is lower by ~ 6 o" and the radius of the convective envelope is larger by ~ 15 a 
with respect to the helioseismic results. Detailed studies have been done to resolve this controversy, see 
e.g.[7|. The latest determinations of the solar photospheric composition Oil] alleviate the discrepancies 
but a definitive solution of the "solar composition problem" still has to be obtained. 

The main effect of changing the heavy element admixture is to modify the opacity profile of the 
sun. It is, thus, evident that the comprehension of the solar composition problem is intimately related 
to understanding the role of opacity in solar modelling. Several authors have investigated the effects 
of opacity changes on the solar structure by using different methods and assumptions (see e.g. [8] and 
references therein). Here, we continue their work, completing and extending the analysis of [H] by using 
a different and original approach. The final goal is to provide the instruments to analyze in transparent 
and efficient way the role of the opacity in the sun and to perform a critical "step-by-step" discussion 
of the present constraints on opacity (and composition) provided by observable properties of the sun. 

In order to calculate the effects of arbitrary opacity changes on the sun, we use the linear solar model 
(LSM) approach, presented in [U and briefly summarized in the appendix. In this approach, the 
structure equations of the present sun are linearized and, by estimating the (variation of) the present 
solar composition from the (variation of) the nuclear reaction rates and elemental diffusion efficiency 
in the present sun, we obtain a linear system of ordinary differential equations that can be easily solved 
and that completely determines the physical and chemical properties of the sun. It was shown in [T] that 
this kind of approach reproduces with good accuracy the results of non-linear evolutionary solar models 
and, thus, can be used to study the role of parameters and assumptions in solar model construction. 

By considering localised opacity changes in LSM approach, we determine numerically the kernels 
that, in a linear approximation, relate an arbitrary opacity variation to the corresponding modification 
of the solar observable properties. These opacity kernels are useful in several respects. First, they 
allow us to individuate the region of the sun whose opacity is probed with maximal sensitivity by each 
observable quantity. Then, they permit us to show that effects produced by variations of opacity in 
different region of the sun can compensate among each others. Finally, they will be used to discuss how 
the different pieces of observational information cooperate to determine the present constraints on the 
opacity profile of the sun. 

The plan of the paper is the following. In sect.[2| we discuss the relation between the effects produced 
by a variation of the heavy element admixture and those produced by a modification of the radiative 
opacity. We show that the relevant quantity is the variation of the opacity profile SK{r) defined in 
eq. Q, that is approximately given by the superposition of the intrinsic opacity change 6Ki{r) and 
the composition opacity change SKz{r), defined in eqs. ([s]) and ([s]) respectively. In sect, [sj we define 
the opacity kernels and we describe the method adopted to calculate them. In sect. |4j we calculate 
numerically the kernels for the squared isothermal sound speed u{r) = P{r) / p[r), the surface helium 
abundance lb, the radius of the convective envelope i?b and the various neutrino fluxes ^i,. Moreover, 
we discuss the constraints on 5K.{r) provided by the present observational data. Finally, in sect, [s] we 
summarize our results. A conclusive view of the constraints on 5K,{r) is provided by fig. [Sl 
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Figure 1: Left Panel: The difference between the squared isothermal sound speed inferred from helioseismic data 
and the predictions of solar models implementing AGS05 (black), GS98 (red) and AGSS09 (blue) heavy element 
admixtures. Right panel: The composition opacity changes SKz{r) that corresponds to using the GS98 (red) and 
the AGSS09 (blue) composition in place of the AGS 05 composition. 



2 The relation between opacity and metals 

We consider a modification of the opacity k{p,T,Y, Z^) and/or of the heavy element photospheric ad- 
mixture {^i}, expressed here in terms of the quantities z\ = Z\^IX\, where Z-^^ is the surface abundance 
of the z— element and X\, is that of hydrogen. If we neglect the role of metals in the equation of state 
and in the energy generation coefficient, the only effect of these changes is to modify radiative energy 
transport in the sun. The relevant parameter, in this respect, is the total variation of the opacity in 
the shell r of the present sun, given by: 

5r}^Hr\ - <Pir),nr),Y{r),Z,ir)) 

^ ^ K{pir),T{r),Y{r),Z,{r)) ^ ^ 

where K{p{r),T{r),Y{r), Zi{r)) is the opacity profile of the SSM, while K{p{r),T{r),Y{r), Zi{r)) is the 
opacity profile of the solar model that implements the modified opacity and photospheric compositior[^ 
We note that the quantity (r) is not related in a direct way to the performed variations of opacity 
and composition. In order to calculate it, we have to take into account that the "perturbed" sun has 
different density (p), temperature (T) and chemical composition profiles {Y and Zi) with respect to the 
SSM. These are not known a priori but have to be obtained as a result of numerical solar modelling. 

A relevant simplification is obtained in the LSM approach presented in ref.[T], where one assumes 
that: i) the performed changes of opacity and heavy element admixture are small; ii) the (variation of) 
the chemical composition of the sun can be estimated from the (variation of) nuclear reaction rate and 
diffusion efficiency of the present sun; Hi) the variation of the metal admixture has a negligible direct 
effect on nuclear production of helium and on diffusion efficiency. In this case, the relation between 
6K^°^{r) and the performed variation of opacity and admixture can be worked out explicitly, as it is 
described in the appendix. It can be shown, moreover, that the source term dK{r) that drives the 
modification of the solar propertied and that can be constrained by observational data can be written 

^ The notation Q indicates, here and in the following, the SSM's value for the generic quantity Q. 

^In mathematical terms, the quantity SK{r) represents the inhomogeneous term in the linearized structure equation of 
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GS98 - AGS05 
AGSS09 - AGS05 


-0.019 
-0.0056 


0.019 
0.0028 


-0.011 
-0.0014 


0.14 
0.018 


0.26 
0.034 


0.56 
0.13 


0.49 
0.12 



Table 1: The variations of helioseismic and solar neutrino observables produced by a variation of the heavy 
element admixture. The above results have been estimated within the LSM approach, by applying the opacity 
changes SKz{r) shown in the right panel of fig^ For CNO neutrinos, we also considered that the fluxes scales 
proportionally to the total CN-abundance. Note that the absolute variation is reported for the surface helium 
abundance, whereas the relative variations are shown for all the other quantities. 



as the sum of two contributions: 

5K{r) = 5Ki{r) + 6Kz{r) (2) 

The first term SKi{r), which we refer to as intrinsic opacity change, represents the fractional variation 
of the opacity along the SSM profile and it is given by: 

^ K(p(r),r(r),F(r),g.(r)) _ 
' Mp{r),T{r),Y{r),Z,{r)) ^ ^ 

This contribution is obtained when we revise the opacity function k(p,T,Y, Zi) and/or we introduce 
new effects, like e.g. the accumulation of few GeVs WIMPs in the solar core (see e.g.|10j and references 
therein) that mimics a decrease of the opacity at the solar center. 

The second term dnzif), which we refer to as composition opacity change, describes the effects of 
a variation of {zi}. It takes into account that a modification of the photospheric admixture implies a 
different distribution of metals inside the sun and, thus, a different opacity profile, even if the function 
k(p,T, Y, Zi) is unchanged. The contribution 5Kz{r) is given by (see appendix A and B): 

^ mr),nr),nr),Z.ir)) _ 
' K{p{r),T{r),Y{r),Z,{r)) 

where Zi(r) = Zi{r) (z\l^\) and can be calculated as: 

51nK 



d\nZ, 



5^i,b (5) 
SSM 



where 5zi represents the fractional variation of Zi and the symbol |ssM indicates that we calculate the 
derivatives along the density, temperature and chemical composition profiles predicted by the SSM. 

Eq.Q, although being approximate, is quite useful because it makes explicit the connection (and 
the degeneracy) between the effects produced by a modification of the radiative opacity and of those 
produced by a modification of the heavy element admixture. In this paper, we take as a reference the 
Asplund, Grevesse, Sauval 2005 composition (AGS05) [2] and we refer with 'SSM predictions' to the 
numerical results obtained by using this composition as input for evolutionary solar model calculation^ 



the present sun. 

^The black dashed lines in the left panel of Fig.l has been obtained by using the FRANEC code. See jllj for a 
description of the code and [T] for a description of the results. 
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Other compilation can be considered, like e.g. the Grevesse, Sauval 1998 (GS98) [5j or the more recent 
Asplund, Grevesse, Sauval, Scott 2009 (AGSS09) [3j. The red and blue dashed lines in the right panel 
of Fig |T] correspond to the opacity change 5K,z{r) that are obtained when we use the GS98 and the 
AGSS09 admixture in place of the AS05 composition, as calculated by applying rel.([5]) and by using 
the logarithmic derivatives dlnK,/ d\\iZ\ presented in Fig. 12 of 0. We observe that the variation from 
AGS05 to GS98 (AGSS09) heavy element admixture corresponds to increasing the opacity by about 
5% (1%) at the center of the sun and by about 20% (5%) at the bottom of the convective region. The 
effects of these opacity changes on helioseismic and solar neutrino observables, calculated in the LSM 
approach, are described in tab{T] and in the left panel of fig{l} They can be compared with the results 
of full non-linear evolutionary codes reported e.g. in [T7j, obtaining a satisfactory agreement. 



3 The method 

In the following, we consider the effects produced by a generic variation of the opacity profile (5K(r), 
without discussing whether this is due to a change of the function T, y, Z^) or to a change of the 
admixture {^^i}. As a results of this modification, we obtain a solar model that deviates from SSM 
predictions. If the opacity variation is sufficiently small (i.e. 5K{r) <^ 1), the sun responds linearly. In 
this case, the fractional variation of a generic quantity Q, defined as: 

«=g-l (6) 

can be related to 5«;(r) by the linear relation: 

5Q = JdrKgir) 5K{r) (7) 

The kernel KQ{r) represents the functional derivative with respect to opacity and allows to quantify 
the sensitivity of Q to opacity variations in different zones of the sun. 

In this paper, we determine numerically the kernels KQ{r) for helioseismic observables and solar 
neutrino fluxes, by using the linear solar model (LSM) approach presented in [H [9] and briefly sum- 
marized in the appendix. In this approach, the structure equations of the present sun are linearized 
and, by estimating the (variation of) the present solar composition from the (variation of) the nuclear 
reaction rates and elemental diffusion efficiency in the present sun, we obtain a linear system of ordinary 
differential equations that can be easily solved and that completely determines the physical and chemi- 
cal properties of the sun. It was shown in [T] that this kind of approach reproduces with good accuracy 
the results of non-linear evolutionary solar models and, thus, can be used to study in an efficient and 
transparent way the role of parameters and assumptions in solar model constructiorj^ 

The estimate of Kq^r) at a given point r = tq is obtained by performing a localized increase of 
opacity in the vicinity of tq. More precisely, we calculate the variation SQ{ro) produced by a normalizecj^ 

*The results presented in this paper can be compared with those presented in |8j where some of the kernels presented 
here are calculated by using static solar models and/or evolutionary models with simplified equation of state and without 
elemental diffusion. Where comparison is possible, a very good agreement is achieved showing that the linearization 
procedure adopted here and the simplifying assumptions implied in fSl do not introduce relevant errors. In order to make 
the comparison, one should note that the definition of the kernels given in 8 differs form that adopted in this paper. 

^We remark that the LSM are linear "by construction". The validity of the the linear relation ([t]) is, thus, not limited 
by the condition SK{r) <C 1. 
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Figure 2: The functions 5Ki^{r), Snouti'i'): Skq^v) and SKi{r) adopted in the parametrizations (11\12 ) and defined 



in eqs.[13) 



gaussian increase of opacity centered in ro: 

5K{r) = G{r — ro) 
with 5r = O.Oli?0, and we assume that: 



exp 



(8) 



(9) 



KQ{ro) = 6Q{ro) 
This corresponds to the approximation: 

Kgiro) ^ j dr KQ{r) G{r - tq) (10) 

which is adequate to describe ah the situations in which we consider opacity variations 5K,{r) which 
vary on scale larger than 6r = O.Oli?©. 

The kernels Kgi^r) can be used to calculate the effects of an arbitrary opacity change 5K{r) and 
allow us to discuss the role of opacity in general terms. In order to consider specific situations and to 
understand what kind of experimental constraints are provided by each observable quantity Q, it is 
useful, however, to consider the simple parametrizations: 



An Sl^inix) + ^out SHontix) 

Aq 6ko{x) + Ai 6ki{x) 



5n{x) 

with A[^, Ao^^t, Aq and Ai free adjustable parameters and: 

Si^outir) = 1 



(11) 
(12) 



1 



exp [(r - rc)/A] + 1 



5Ko(r) = 1 



6Ki{r) 



Rh 



(13) 
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dQ 




SQout 


6Qo 


6Qi 


6u{0.1Rq) 


0.019 


-0.036 


-0.017 


-0.022 


5u{0.2Rq) 


0.052 


-0.054 


-0.0025 


-0.014 


Su{OARq) 


-0.084 


0.087 


0.0037 


0.031 


6u{0.65Rq) 


-0.16 


0.17 


0.011 


0.11 


AFb 


0.073 


0.069 


-F0.142 


0.062 


5Rb 


0.12 


-0.14 


-0.02 


-0.10 




-0.069 


-0.031 


-0.100 


-0.030 


S^Be 


0.85 


0.41 


1.26 


0.38 




1.93 


0.75 


2.68 


0.68 




1.65 


0.50 


2.15 


0.48 




1.14 


0.28 


1.43 


0.30 



Table 2: The coefficients SQm, SQout, SQq and SQi defined in eq.(H), which allow to calculate the response of 



helioseismic and solar neutrino observables to opacity changes parameterized by eqs. 



where Tc = O.3i?0, A = 0.01, while Rh = 0.7307?© is the radius of the convective envelope predicted by 
SSM (see [1|). The functions SKiri{r) and SKo^tir) correspond to a constant increase of the opacity in 
the energy producing zone (r < 0.3Rq) and in the outer radiative region (r > O.3i?0), see Figj2j They 
have been defined in such a way that S Kin{r) + 6 Hout{r) = 1- The function 6Ko{r) and 5Ki{r) correspond 
to a global rescaling and to a linear tilt of the opacity profile. 

In linear approximation, the fractional variation 6Q produced by the opacity profiles ( 11|12 ) can be 
expressed as: 



5Q = 
SQ = 

where the coefficients 6Qj are given by: 



Ain 5Qin -\- j4out <5Qout 

Aq 5Qo + Ai 6Qi 



dr Kg^r) 6Kj{r) 



(14) 



(15) 



with j = in, out, 0, 1. We report these coefficients in Tabj2] for helioseismic observables and solar 
neutrino fluxes. 



4 The opacity kernels 

We calculate the opacity kernels for following observable quantities: the squared isothermal sound 
speed u{r) = P{r)/p{r); the surface helium abundance lb; the depth R\, of the convective envelope; the 
solar neutrino fluxes ^i^, where the index v = pp. Be, B, N, O refers to the neutrino producing reactions 
according to the usual convention. Our results are presented in the following sub-sections, together 
with a discussion of the present observational constraints on 6Hi{r). 



7 




4.1 The sound speed 

In Figjsj we discuss the properties of the sound speed kernel Ku{r, r') defined by: 

6u{r) = / dr' Ku{r,r') 5K{r') (16) 



where 5u{r) is the fractional variation of the squared isothermal sound speed u{r) = P{r) / p{r). In 
the left panel of Fig. [sj we show the functions fr'{r) = Ku{r, r') Rq, calculated for the selected values 
r' = 0.1, 0.2, . . . , O.7i?0. In our approach, the functions fr'ir) correspond to the sound speed variations 
produced (in a linear theory) by the localised opacity increases G{r — r'). One has large effects close to r' 
which are due to the variation of the temperature profile of the sun, as can be understood by considering 



that 5u = dP — dp ~ 6T — Sp., where fj, is the mean molecular weight of the solar plasma. From eq. (57), 
we see that 6T[r) is expected to have a sharp decrease close to r', by an amount approximately equal 
to l//((r'), while 5/i remains approximately constant, being related to the chemical composition of the 
solar plasma. 

In the right panel of Fig. [s] , we show the functions grir') = Ku{r, r') Rq, calculated for the selected 
values r = 0.1, . . . , 0.7 Rq, that quantify the sensitivity of the sound speed at a given r to the opacity in 
the shell r' of the sun. They clearly indicates that u{r) is maximally sensitive to the value of the opacity 
at r' ~ r. However, the displayed functions are different from zero everywhere, and have negative values 
in large part of the solar radiative region. This suggests that compensating effects can occur, especially 
when one considers opacity modifications that extend over a broad region of the sun. 

In this respect, it is important to note that the sound speed is practically insensitive to a global 
rescaling of opacity. This can be appreciated by looking at the black dotted line in the left panel of 
Fig. |3] , which is defined by: 

6uo{r) = dr'Ku{r,r') (17) 



and corresponds to the sound speed variation produced (in a linear theory) by a constant rescaling 
SKo{r) = 1. For visualization purposes, it is useful to note that the values of 6uo{r) at r = 0.1, . . . , 0.7Rq 
correspond to the integral in r' /Rq of the functions gr{r') displayed in the right panel of Fig. [sj We see 
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Figure 4: Left Panel: The functions Suinir), Suoutif), ^^Q^r) and 5ui{r) defined in eq.(18), which allow to 

ference 



calculate the sound speed response to opacity changes parametrized by eqs.{ ll\T2 ). Right panel: The diffe 
between the sound speed inferred from helioseismic data and that obtained by solar models implementing AGS05 
heavy element admixture (dotted line) with suitably modified opacity profile (the solid lines correspond to models 
A, B, C, D described in the text). 



that 5uQ{r) is very small, as a result of an almost perfect compensation between the positive contribution 
from the region r' ~ r and the negative contribution from the other regions of the sun. 

A qualitative argument to explain the stability of the sound speed is the following. The virial 
theorem connects the gravitational energy, Eg = — J dm Gm/r, and the thermal energy content of a 
given star, Ei = 3/2 / dm {P/ p) = 3/2 j dm u, being Eg = —2Ei. For a generic star, a global rescaling 
of the opacity reflects into a global rescaling of the radial profile and, thus, taking into account the 
virial theorem, also into a global rescaling of the sound speed. This is not the case for the sun because 
the solar radius is observationally determined. We are forced to re-adjust the free parameters in the 
model in order to keep the solar radius fixed, with the effect of stabilizing the radial profile and, thus, 
the sound speed profile u{r). 

The above result has relevant implications. In particular, the statement that the "sound speed 
problem" requires an increase of the opacity at the bottom of the convective region is, strictly speaking, 
not correct because other solutions are possible. We can consider, e.g, the sound speed profiles produced 
by the opacity variations parametrized by eqs. ( 11|12 ). In linear approximation, we have that: 



6u{r) = Ai^duinir) + AontSuontir) 

5u(r) = Aq duo{r) + Ai 5ui{r) . (18) 

The functions 6uia{r), Suoutir), Suo{r) and 6ui{r) are shown in the left panel of Fig. [4] and are reported 
in Tab. [2] for the selected values r = 0.1, 0.2, 0.4, 0.65 i?©. We see that Suo{r) ^ 6ui{r), indicating 
that the "tilt" and not the scale of the opacity is fixed by the sound speed. Moreover, we have 
6uin{r) — —5uout{r) that shows that the effect of an enhancement of the opacity in the external 
radiative region can be equally produced by a decrease of the opacity at the solar center. 

This is confirmed in the right panel of Fig|4] where we show the sound speed profiles obtained by 
implementing four different opacity modifications, described by: 

Model A: 15% decrease of opacity in the energy producing region (^in = —0.15 and ^out = 0); 
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Model B: 15% increase of the opacity in the outer radiative region {A\a = and ^out = 0.15); 

Model C: a hnear tilt of opacity corresponding to a 20% increase at the bottom of the convective 

envelope (^o = and Ai = 0.2); 

Model D: a linear tilt plus a global rescaling of opacity corresponding to a 20% decrease at the solar 
center (^o = -0.2 and Ai = 0.2). 

The sound speed profiles obtained in all the considered cases reproduce equally well the result inferred 
by helioseismic data, showing that the helioseismic determination of u{r) translates into a bound on 
the differential increase ^out — An — 0.15 and on the tilt Ai ~ 0.2, with no relevant constraint on Aq 
(or, equivalently, A\^ + Aout) that fix the global scale of opacity. 

In light of this observation, it is intriguing the possibility that non-standard effects that mimic a 
decrease of the opacity at the solar center, like e.g. the accumulation of few GeVs WIMPs in the solar 
core (see e.g. |10]). could have a role in the solution of the solar composition puzzle. We will see, in 
the following, that this possibility is disfavoured by the determination of the surface helium abundance 
and by the measurement of the ^Be and neutrino fluxea^ 



4.2 The surface helium abundance 

In the left panel of Fig jsj we show with solid line the functional derivative Ky (r) of the surface helium 
abundance Y\^, defined according to equation: 

AFb = J dr Kyir) 6K{r) (19) 

where we considered the absolute variation AYb to conform with the notations adopted in [iQ 

The surface helium abundance depends on the initial chemical composition and on the effects of 
elemental diffusion. In ref.^, by assuming that (the variation of) the time-integrated effect of diffusion 
can be estimated from (the variation of) diffusion efficiency in the present sun, we have obtained the 
following relation: 

AFb = Ay AFini + Ac SC (20) 

with Ay = 0.838 and Aq = 0.033, which gives AYb as a function of the absolute variation of the initial 
helium abundance, AYini, and of the fractional variation of pressure at the bottom of the convective 
region, 6C = SP\,. The above equation allows us to obtain the functional derivative Ky{r) as the sum 
of two contributions, one related to the term Ay Alini and the other to Aq 6C. These are shown in 
FigjS] with red and blue dashed lines, respectively. 

The function Kyir) is positive everywhere, showing that an increase of opacity in an arbitrary 
shell of the sun translates into an increase of the helium abundanc^ Moreover, the kernel Ky{r) 
has a rather broad profile. This implies that the determination of lb effectively constrains the opacity 
scale and breaks the degeneracy between the possible solutions of the sound speed problem presented 



^The idea that the accumulation of few GeV WIMPs in the solar core could alleviate the "solar composition problem" 
was originally proposed by |13| . The recent paper |14| presented a qualitative implementation of this idea. The effect on 
boron neutrinos is discussed in |15| . 

^Here and in the following, we use the notation Qh to indicate that a given quantity Q{r) is evaluated at the bottom 
of the convective region, i.e. Qh = QiRh) where _Rb ~ 0.730 7?©. 

* If the chemical composition is fixed, an increase of the opacity implies a decrease of the total luminosity L which 
roughly scales as L oc (J-^/k, where fi is the mean molecular weight. In order to reproduce the observed solar luminosity 
Lq, we are forced to readjust the chemical composition of the sun by increasing the helium abundance. This has the 
simultaneous effects of increasing ^ and decreasing k. 
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Figure 5: Left Panel: The solid line corresponds to the kernel Ky^r) defined in eq.(19). The red and blue 
dashed lines describe the contributions to surface helium variations provided by the terms ylyAi^n; and AqSC, 
respectively. See text for details. Right Panel: The solid line corresponds to the kernel Kji(r) defined in eq.(23). 
The red, blue and green dashed lines describe the contributions to convective radius variations provided by the 
terms FyAl^ni, TcSC andTi^dK^,, respectively. See text for details. 



in the previous section. The SSM that implements the AGS05 heavy element admixture predicts 
the value Y]^ = 0.229 (see e.g. |lj) which is about 6a lower than the helioseismic determination 
lb = 0.2485 lb 0.0034 [7]. This discrepancy requires an increase of the opacity, as it can be obtained, 
e.g., by increasing the metal content of the sun. Models that accounts for a reduction "effective" opacity 
at the solar center are expected to decrease Yb, increasing the disagreement with helioseismic results. 

A simple quantitative analysis can be performed by considering the opacity profiles ( 11|12 ) that 
produce the variations Al^ given by: 

AFb = 0.073 An + 0.069 Aout =^ 0.07 (An + ^out) (21) 
AFb = 0.142^ + 0.062 ^1 (22) 

By comparing the SSM prediction with the helioseismic result, we obtain + Aut = 0.28 it 0.05 
that, combined with the "orthogonal" constraint ^out — ^in — 0.15 provided by the sound speed, gives 



Ain = 0.04 H- 0.09 and ^out = 0.19 -r- 0.24. Alternatively, we can use eq.(22) and the information 
on the tilt Ai ~ 0.2 provided by the sound speed, to fix the opacity at the solar center obtaining 
Aq = 0.035^0.075. 



4.3 The convective radius 

In the right panel of Figjsj we show with a solid line the functional derivative Kfi(r) of the convective 
radius -Rb defined according to the relation: 

5i?b = Jdr Kn{r) 6K{r) (23) 

We see that the kernel Kji{r) has a very sharp peak at r ~ i?b = O.73O-R0 that reflects the (well-known) 
fact that the convective radius is particularly sensitive to the opacity at the bottom of the convective 
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region. The shape of this peak has not a precise physical meaning and depends on the method of 
calculation. The effect of opacity changes has been, in fact, estimated within the LSM approach, in 
which the fractional variation 5Rh is calculated from: 



(24) 



where Fy = 0.449, Tc = —0.117 and = —0.085, while (5Kb is fractional variation of opacity at 
the bottom of the convective envelope, i.e. 5k]^ = 6k{R]^). The "local" term r^^Kb translates into a 
delta-function contribution 6{r — i?b) to the functional derivative. Since we evaluate numerically 
the kernel by applying a localised gaussian increase of opacity, this is convolved with the function 
G{r — ro) given in eq.([8]). As a final result, one obtains the contribution F^ G{r — R\^) to the kernel 
Kji[r) which is shown by the green dashed line in Figjsj whereas the red and blue dashed lines describe 
the contributions arising from Fy Alini and Tc SC, respectively. We remark that the area under the 
peak at r = i?b is approximately equal to F^ and does not depend on the calculation method. We can, 
thus, safely use the functional derivative Kfi(r) to describe all the situations in which opacity varies on 
scale larger than 6r = O.Oli?0. 



Eqs.(20) and (24) can be combined to obtain a direct determination of (5Kb from quantities that are 



all determined by helioseismic observations. We can, in fact, eliminate Alini from eq.(24), obtaining: 



(5Kb = Cy AFb + Cr 6Rb + Cp 5p^ 



where: 



Cy 
Cr 

Cn 



Ty 

Ay F« 



6.27 



1 
1 



= -11.71 

Ac Ty 
Ay 



c 



-1.58 



(25) 

(26) 
(27) 
(28) 



In the derivation of the above result, we have considered that 5C = (5Pb — ^Ph^ since the fractional 
variation of the sound speed 6u{r) is expected to vanish at the bottom of the convective region, i.e. 
(5iib = (5Pb — ^Ph — 0, as it is discussed in [T]. The discrepancy between the helioseismic determinations 
of -Rb and Yb and the predictions of SSMs implementing AS05 admixture is quantified as 5Rh = 
-0.0205 ± 0.0015 and A^b = 0.0195 ± 0.0034. The density at the bottom of the convective region 
deviates from value inferred by helioseismology by 6ph = 0.08, as it is discussed e.g. in [12]. We obtain 
(5Kb — 0.24 lb 0.03, where errors have been combined in quadrature and we have neglected the (sub- 
dominant) contribution to the total error budget due to the uncertainties in the density determination. 
We remark that the obtained results is model- independent, since it does not rely on any assumption or 
parametrization for the function 5K{r). 

As a final application, we consider the response of i?b to the opacity changes parametrized by 
eqs.( 11|12 ). We obtain the relations 



6Ri, = 0.12 An -0.14 Aout 
SR^ = -0.02 Aq - 0.10 Ai 



0.13 (An - A 



out ) 



(29) 
(30) 



which show that the convective radius, just like the sound speed, provides bounds on the differential 
increase j4out — An and on tilt Ai, with no relevant constraints on Aq (or equivalently + Aut) that 
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Figure 6: Left Panel: The solar neutrino kernels K^{r) defined in eq.(31). Right Panel: The solid lines are 



the normalized solar neutrino kernels Ki,{r) / S'^nfi. The dashed line shows the normalized kernel KT{r) / ST^fi 
defined in eq. \3^, that describes the response of the solar central temperature to localised opacity modifications. 



fix the global scale of opacity. By considering dR^ = —0.0205 it 0.0015, we obtain ^out — ^in ~ 0.15 
and Ai ~ 0.2 in substantial agreement (and complete degeneracy) with the information provided 
by the sound speed measurement. By combining this information with the constraints provided by 
the surface helium abundance and performing a simple analysis, we obtain Ai^ = 0.07 it 0.04 and 
^out = 0.21±0.04 for the parametrization given in eq. (fTT| and = 0.056±0.040 and Ai = 0.187±0.023 



for that given in eq. (12). The corresponding bounds on the opacity change 6K{r) are shown in fig. 
and commented in the conclusive section. 



4.4 Neutrino Fluxes 

In the left panel of Figj6| we show the functional derivatives Ky{r) of the neutrino fluxes defined 
according to relation: 

= j dr Ku{r) 6K{r) (31) 

where the index u = pp, Be, B, N, O labels the neutrino producing reactions according to the usual 
convention. The kernel Ku{r) have been calculated by using the LSM approach and by taking into 
account that the fractional variations of the fluxes are related to physical and chemical properties 
of the sun by: 

5^, = Jdr [0,,p(r) 6p{r) + </),,T(r) 5T{r) + </.,,y(r) AY{r) + 0,,z(r) 6Z{r)] (32) 

The functions (pujif) have been defined and calculated in the fig. 10 of [1]. 

Our results show that neutrino fluxes probe the opacity of the sun in the region r < O.45i?0. The 
kernels Kb{x), K^d^), Kt<}{x) and Ko{x) are positive- valued almost everywhere, while the kernel 
Kpp{r) is negative. This indicates that an increase of opacity generally translates into an enhancement 
of '^B, ^Be and CNO neutrino fluxes and into a (slight) decrease of the pp— neutrino component. In 
Tab. [2| we show the coefficients 6^u,in, 5^u,out, 5^u,o and 6^^,! that allows to describe the effects of 
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opacity changes parametrized by eqs. ( 11|12 ), through the simple relations: 



6<P 



V 



^in ^^v^r, + ^out ^^v 



(33) 
(34) 



We see that and CNO neutrinos are extremely sensitive to opacity changes, as it is expected since 
they strongly depend on the temperature of the central regions of the sun. It is interesting to note in 
the right panel of Fig. [ojthat the normalized neutrino kernels Ky{x) j^yi^ have a common behaviour, 
with two maxima at r ~ O.l-R© and r ~ O.3-R0, and one minimum at r ~ 0.2i?Q. This shows that 
different fluxes basically probe the same quantity, constraining the opacity profile with the maximal 
sensitivity in the two regions of the sun r ~ 0.05 0.15 and r ~ 0.2 0.45. 

Some insights on the above results can be obtained by recalling that the total neutrino flux is 
essentially fixed by the solar luminosity constraint. We, thus, expect that A<^tot = Yliv '^v ' — 
0, where are the SSM predictions for the various neutrino components. Considering that about 
99% of the total flux is provided by pp and Be neutrinos, the luminosity constraint implies <5$pp ~ 
— ('^'Be/^pp) <5^Be = — 0.075 (5^>Be which explains the smallness of the pp-neutrino kernel, the ratio 
between the pp and Be-neutrino coefficients in tab. [2] and the equality between the normalized pp and 
Be-neutrino kernels observed in the right panel of Fig. |6j The wavy shape of the kernels Ky{r^ reflects, 
instead, the response of central temperature to localised opacity modifications. This was first noted 
and discussed by |8] and it is seen in the right panel of Fig. [6j where we show with a dashed line the 
functional derivative of the central temperature Tc, defined by: 



For convenience, we plot the normalized kernel K'j^ix^jbTci^^ where the normalization factor is ^T^^ = 



The peculiar behaviour with r of the kernel Kxir) cannot be explained in simple terms. A qualitative 
comprehension can be obtained from Fig. [7| where we show the temperature profiles 5T{r) produced by 
the opacity changes 5K,{r) = G{r — vq) with tq = 0.1, 0.2, 0.3 Rq. We see that the performed opacity 
modifications translate into a large increase of temperature close to ro, that necessarily alters nuclear 
burning rates since tq is inside or close to the energy producing region. The free parameters of the solar 
model are readjusted in such a way that the same luminosity is obtained with a different temperature 
profile. 

In the two cases ro = O.IRq (i.e. opacity increase well inside the energy producing region) and 
ro = O.3i?0 (i.e. opacity increase just outside the energy producing region), the entire energy producing 
zone is affected. A new "equilibrium" situation is achieved in which nuclear burnings occur at higher 
temperatures with a larger helium abundance, favoring ^B, CNO and ^Be neutrinos at expenses of 
pp-neutrinos. For ro = 0.2Rq, we have a peculiar situation since the maximal effect is produced in a 
region of the sun where only pp neutrinos are produced, as it is seen from fig. 10 of |1]. This necessarily 
shifts energy production outwards with respect to the SSM case. In order to avoid overproduction of 
energy and to recover the observed luminosity, the central temperature Tc has to be (slightly) decreased 
so that the burning rates at the center of the sun are suppressed. 

At present, the best studied components of the solar neutrino flux are the ^B neutrino flux which 
is determined by the SNO neutral current measurement with about 6% accuracy, <I>b = (5.18 it 0.29) x 
10^ cm~^ s~^ |16| . and the ^Be neutrino flux which is measured by Borexino, ^>Be = (5.18 it 0.51) x 
10^ cm~^ s~^ [T7]. These fluxes have to be compared with the results of theoretical calculations. SSMs 




(35) 



0.138. 
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implementing AGS05 heavy elements admixture predict values for $b and <I>Be which are about 10% 
lower than the experimental results. We take, for definiteness, the values $b = 4.66 x 10^ cm~^ s~^ 
and $Be = 4.54 x 10^ cm~^ s~^ obtained in jT2] which are affected by ~ 9% and ~ 5% theoretical 
uncertaintie^ respectively. The difference between the theoretical predictions and the experimental 
data points toward a moderate increase of the central opacity of the sun. As an example, a 5% 
increase of the opacity in the region r < 0.3Rq would produce a 9.7% (4.3%) increase of the *B (''Be) 
neutrino flux. Models that account for a reduction of the "effective" central opacity, due e.g. to WIMP 
accumulation, increase the disagreement and are, thus, disfavoured by solar neutrino data. 

We, finally, discuss the constraints provided by solar neutrinos on opacity changes parametrized by 
Eqs. ( 11|12 ). We note that the coefficients Am and ^out and/or Aq and Ai required to fit helioseismic 
results (see the previous section), produce enhancements of the and ^Be neutrino fluxes which 
are equal to 28% and 14%, respectively. While the ^Be component would be consistent with the 
observational data, the ^B neutrino flux is too large with respect to the SNO measurement. When 
we fit simultaneously helioseismic and solar neutrino data, we obtain a slight reduction of the required 
opacity change. A simple x^-analysis gives = 0.05±0.03, ^out = 0.19±0.03 for the parameterization 
( [IT| ) and Aq = 0.038 ± 0.034, Ai = 0.192 ± 0.023 for the parameterization ^ with the best fit values 
corresponding to Xmin/'i-°-f- — 2.1/2 and Xmin/^-°-f- — l-V^, respectively. The corresponding bounds 
on 6K{r) are displayed in fig. Isl 



5 Summary 

In this paper, motivated by the solar composition problem and by using the recently developed LSM 
approach p], we have discussed the effects of arbitrary opacity changes on the sun. Our main results 
can be summarized as it follow: 

i) We have discussed the relation between the effects produced by a variation of the heavy element 
admixture and those produced by a modification of the radiative opacity. We have shown that the 

^The quoted theoretical uncertainties are obtained from [THj. In our estimate, we have not included the contribution 
due to opacity, since this is considered as a free parameter in our analysis. 
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Figure 8: Left Panel: The constraints on the opacity profile SK{r) obtained from helioseismic data. Right panel: 
The constraints on the opacity profile Sn^r) obtained from helioseismic and solar neutrino data. See text for 
details. 



relevant quantity is the variation of the opacity profile 6K{r) defined in eq.Q, that is approximately 
given by the superposition of the intrinsic opacity change 6Ki{r) and the composition opacity change 
6Kz{r), defined in eqs.Q and ^ respectively. 

ii) We have studied the response of the sun to an arbitrary modification of the opacity 5K{r). Namely, 
we have calculated numerically the kernels that, in a linear approximation, relate the opacity change 
5K{r) to the corresponding modifications of the solar observable properties. We have considered the 
following observable quantities: the squared isothermal sound speed u{r); the surface helium abundance 
Yb; the convective radius R\^; the solar neutrino fiuxes (^y. 

Hi) We have shown that different observable quantities probe different regions of the sun. Moreover, 
effects produced by variations of opacity in distinct zones of the sun may compensate among each other. 
In this respect, we noted that the sound speed u{r) and the depth of the convective envelope are 
practically insensitive to a global rescaling of the opacity. 

iv) As a consequence of the above result, we have seen that the discrepancy between the SSM predictions 
for u{r) and i?b and the helioseismic inferred values can be equally solved by a ~ 15% decrease of the 
opacity a the center of the sun or by a ~ 15% increase of the opacity in the external radiative region. The 
degeneracy between these two possible solutions is broken by the "orthogonal" information provided 
by the measurements of the surface helium abundance and of the boron and beryllium neutrino fluxes, 
that fix the scale of opacity and indicate that only the second possibility can effectively solve the solar 
composition problem. 

v) We have derived a model-independent relation, eq. (25), that allows to obtain a direct determination 
of the opacity at the bottom of the convective envelope 5k\, from quantities that are determined by 
helioseimic observations, i.e. the surface helium abundance AYi3, the convective radius 5R\j and the 
density at the bottom of the convective region dp^,. By considering the present observation values, we 
have obtained 5K\y = 0.24 it 0.03. 

A conclusive view of the present constraints on 5K{r) is contained in Fig. [sj The red (blue) dashed 
lines correspond to the composition opacity changes obtained when we replace the AGS05 composition 
with the GS98 (AGSS09) admixture. The blue tick at r = i?b is the value of the opacity at the bottom of 
the convective region that is obtained by applying the model-independent rel.(25) to helioseismic data. 
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The dark and light areas individuate the opacity changes 6K(r) that, for each value of r, are obtained 
at 68.3% (Icr) and 95.4% (2a) C.L. by applying a simple analysis to the opacity modification 
parametrized by eq.(12). The left panel takes into account only helioseismic observables. Namely, it 
is obtained by using the present observational determinations of the surface helium abundance and 
of the convective radius. The selected opacity profile 5K{r) is also consistent with the helioseismic 
determination of the sound speed, as it is has been discussed in fig. |4j The right panel includes also 
the observational information on boron and beryllium neutrinos. These data move the required opacity 
change 6K(r) towards slightly smaller values. In summary, we see that the helioseismic and solar 
neutrino data select an opacity profile that well corresponds to the GS98 heavy element admixture. 
They clearly disfavour different solutions of the solar composition problem, like e.g. that recently 
proposed by [14j . 
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Appendix A: The total variation of opacity (5/ttot(^) 



In order to calculate the total variation of opacity SK,tot{r) at a given radius r, we have to take into 
account that the perturbed sun has different temperature, density and chemical composition profiles 
with respect to SSM. We define: 



K{p{r),T{r),Y{r),Zi{r)) 



K{p{r),T{r),Y{r),Z,{r)) 
and, by expanding the above equation to linear order, we obtain: 

5K}°\r) = KT ST{r) + Kp 5p{r) + ky ^Y{r) + ^ 5Z;{r) + 5ni{r) 



(36) 



(37) 



where: 



KY{r) 



91nK 



d\nT 
(91n K 
dY 



SSM 



SSM 



Ki(r) 



d\np 
dlnZi 



SSM 



SSM 



and the symbol |ssM indicates that we that we calculate the derivatives Hj{r) along the density, tem- 
perature and chemical composition profiles predicted by the SSM. The quantity 5Ki{r) represents the 
intrinsic opacity change and it is given by: 



5Ki{r) 



K{p{r),Tir),Y{r),Zi{r)) 



1 



(38) 



K{p{r),T{r),Y{r),Zi{r)) 

By taking advantage of the equation of state of the stellar plasma, we can eliminate the density from 



Eq. (37). We use the relation: 



where i-V(r) 
obtain: 



Sp{r) = 6P{r) - 5T{r) - PYir)AY{r) (39) 
-5/[8 — 5y(r)] and p represents the mean molecular weight, and we 



—d\n p/dY 1 

5K*°*(r) = (kt - >^p) ST{r) + 6P{r) + {ky - Py «p) Ay(r) + ^ 5Z;{r) + 5Ki{r) (40) 



In order to evaluate eq.(|40[), we need to estimate the chemical composition of the perturbed sun, i.e. the 
quantities 5Z\{r) and /\Y{r). By using the approximate method introduced in ref.fT] and reviewed in 
the next section, we can relate the chemical composition profiles to the modification of the photospheric 
heavy element admixture, to the present values of 5T[r) and 5P{r) and to the parameters AYini and 
6C which represent the absolute variation of the initial helium abundance and the fractional variation 
of the pressure at the bottom of the convective region, respectively. We obtain: 



6K^°\r) = 4 6T{r) + ^ 6P{r) + Ky AYini + kcSC+ [SKi{r) + 6Kz{r)] 
where 5Kz{r) is the composition opacity change given by: 

5Kz{r) = ^Ki5zi 



(41) 



(42) 
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where 5z[ is the fractional variation of Zi = Zi^]^/X]^, while: 

= KT - Kp (1 + Py Ct) - Ky 
Kp = Kp{l- Py Cp) - Ky Cp 
Ky = {ky - Kp Py) + Qy 
Kc = Qc kz 

with the coefficients Qh and defined in the next section. In the derivation of the above equation, 
we took into account that i^i = /^z where kz = dlnK/dlnZ is the partial derivative of opacity with 
respect to the total metal abundance (i.e. calculated by rescaling all the heavy element abundances by 
a constant factor, so that the metal admixture remains fixed). 
It is useful to define the opacity change (5«;(r) given by: 

dK{r) = 6Ki{r) + 6Kz{r) (43) 

which groups together the contributions to (^K*°*(r) which are directly related to the variation of the 
input parameters. While the other terms in eq.(|41[) represent derived quantities that are determined 
by solving the structure equations and/or by fitting the observed properties of the sun, the opacity 
change 6K{r) can be varied, in principle, in arbitrary way. 



Appendix B: The chemical composition of the sun 

The chemical composition of the perturbed sun should be calculated by integrating the the perturbed 
structure and chemical-evolution equations starting from an ad-hoc chemical homogeneous ZAMS 
model. In ref.|l], we proposed a simplified approximate procedure that allows to estimate with sufficient 
accuracy the helium and metal abundances of the modified sun, without requiring to follow explicitly 
its time-evolution. We review this procedure and we extend it to take into account the effect of a 
variation of the photospheric composition. 

In order to quantify the relevance of the different mechanisms determining the present composition 
of the sun, we express the helium and metal abundance according to: 

y(r) = Yini [l + DY{r)]+Y,,^c{r) 

Zi{r) = [l + Dz{r)] (44) 

Here, lini and Z[^i^[ are the initial values for the abundances, the terms DY{r) and Dz{r) describe the 
effects of elemental diffusion and l^uc('") represents the total amount of helium produced in the shell 
r by nuclear processes. We note that we, implicitly, assumed that heavy elements have all the same 
diffusion velocity by introducing a common diffusion term Dz{r) for all metals. 

We are interested in describing how the chemical composition is modified when we perturb the SSM. 
In the radiative core (r < Rh)-, we neglect the effects produced by variations of the diffusion terms 
and we write: 
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Ayfr) = AK, 



1 + Dy(r) +Aynuc(r) 



5Z,{r) = 5Z,,^, (45) 



^"The diffusion terms Dyir) and Dz{r) are at the few per cent level in the radiative region. Their variations are, thus, 
expected to produce very small effects on the solar composition. 
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where AYtiuc(^) is the absolute variation of the amount of hehum produced by nuclear reactions. A 
better accuracy is required in the convective region, because the surface helium abundance Yj, is an 
observable quantity. We, thus, discuss explicitly the role of diffusion and we write: 

5Zi,b = .5Zi,ni + ^ ^ff 5Dz,^ (46) 

t + i^Z.b 

where 5-Dy,b and SDz,b are the fractional variations of the diffusion terms Dy ]^ and -Dz,b- The quantities 
Alb and 6Z]^ are related among each other, since the metals-to-hydrogen ratios at the surface of the 
sun are observationally fixed. If we indicate with zi = Z\^lX\y the surface abundance of the i— element 
(rescaled to that of hydrogen), we obtain: 

5^i,b = - — ^ AFb + (47) 
where we considered that Xb ~ 1 — lb) while bzi is defined by: 

The above relation can be rewritten in terms of the initial helium and metal abundances, obtaining: 

bZi;„^ = Qo AUni + Ql SDY,h + Q2 SDz,^ + Sz; (49) 

The coefficients Qi have been calculated explicitly in ref.[lj and are given by Qq = —1.141, Qi = 0.041 
and Q2 = +0.118, respectively. 

Up to this point, the derived relations have a general validity, since the only assumption implied by 
our analysis is that the heavy elements have all the same diffusion velocity (we take iron as representative 
for all metals). To complete our calculation, we have to estimate the term AYnuc(?') in eq.(45) and 
the quantities SDy^h and SDz,h in eqs.(46, 49). We use the procedure adopted in the LSM approach, 
where we assumed that the helium produced by nuclear reactions scales proportionally to the energy 
generation coefficient (and, thus, the helium production rate) in the present sun, i.e: 

AY^^, = Ynnc{r)6e'^\r) (50) 

where 5e^°^{r) is the fractional variation of the energy generation rate. The effect of elemental diffusion 
is modelled by assuming that the terms D^ ]^ vary proportionally to the efficiency of diffusion in the 
present sun, obtaining (see sect. 6. 3 and appendix C of ref.[lj): 

SDY,b = HyJTb+npJPb 

5Dz,b = Uz6n + Up6Pb (51) 

where Ily = 2.05, Uz = 2.73 and Up = -1.10. 

By following the calculations described in sect. 6.2 of [I], we obtain the following expression for the 
variation of the helium abundance in the radiative region: 

AY{r) = ey (r) Al^^i + CHr) 6T{r) + ^p{r) 6P{r) (52) 

The coefficients ^,hiR) are defined in eq.(31) of ref. [IJ and are shown in their fig. 4. 
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By taking into account relations (51) and by considering the conditions that hold at the bottom of 
the convective region (expressed in eq. (21) of [IJ), we can estimate the variation of metal abundances 
in the radiative region, obtaining: 

5Z,{r) = 5Z,^^, = Qy Allni + Qc6C + 6zi (53) 

where Qy = —0.887, Qc = —0.164 and 5C = 6Ph represents the variation of pressure at the bottom of 
the convective envelope. 

Finally, we can calculate the abundances in the convective region obtaining 

6Y^ = AyAYi^i + AcSC (54) 
(5Zi,b = BYAYir,i + Bc5C + 5zi (55) 

where Ay = 0.838, Ac = 0.033, By = -1.088 and Be = -0.043. 

We remark that, while the quantities AYini and SC are parameters which are univocally determined 
by imposing the appropriate integration conditions (see next section) , the quantities 5zi represent input 
parameters for solar model calculations. 

Appendix C: Linear Solar Models 

By expanding to linear order the structure equations of the present sun close to the SSM solution 
and by assuming that the variation of the chemical abundances of the sun can be estimated by the 
procedure outlined in the previous section, we obtain a linear system of ordinary differential equations 
that completely determine the physical and chemical properties of the "perturbed" sun (see pLj for 
details). Namely, we obtain: 

[jp6P + -fT6T-Sm + jy AYini] (56) 
[(7P -1) 6P + jTST + 6m + -fy AYi^i] 



dSm 


1 


dr 




d6P 


1 


dr 


Vp 


d6l 


1 


dr 


k 


d5T 


1 


dr 


h 



[a'p 6P + aT5T + 5l + a'y AY.^i + a'c SC + 5k\ 

The coefficients 7/j, and and the scale heights 1^ = dln{h)/dr have been calculated in [1] and 
are shown in their fig. 1, fig. 5 and fig. 6. The parameters Al^ni and 6C represent the absolute variation 
of the initial helium abundance and the relative variation of pressure at the bottom of the convective 
envelope and can be univocally determined by imposing the appropriate integration conditions. At the 
center of the sun (r = 0) we have: 

Sm = 7p,o 5Po + 7T,o STq + -fy^ AFini 5P = 5Po 

61 = (3'p^o^Po + (3'Tfl5To + ^'y^o^Yini + Pc,o^C! 6T = 6To 

where the subscript "0" indicates that a given quantity is evaluated at r = 0. At the bottom of the 
convective envelope (r = R\j), we have instead: 

6m = —rficonv 6C 6P = 6C 

51 = 6T = A'y AYini + A'c SC 
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where Ny = 0.626 and A!(. = 0.025 and 

^conv — ^ cam/ / — 0.0192 is the fraction of solar mass 

contained in the convective region. 

The term 5n{r) contains the contributions to the modification of the opacity profile of the sun that 
are directly related to the variation of the input parameters. It is given by: 

5K{r) = 5Ki{r) + 5Kz{r) (57) 

where 5Ki{r) and 6Kz{r) are the intrinsic and composition opacity changes, defined in eqs. ^ and (U 
respectively. It represents the source term that drives the modification of the solar properties an that 
can be bounded by observational data. 
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